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We study the interplay between Coulomb blockade and the Kondo effect in quantum dots. We use 
a self-consistent scheme which describes mesoscopic devices in terms of a collective phase variable 
(slave rotor) and quasiparticle degrees of freedom. In the strong Coulomb blockade regime, we re- 
cover the description of metallic islands in terms of a phase only action. For a dot with well-separated 
levels, our method leads to the Kondo effect. We identify the regime where a crossover between the 
Coulomb blockade regime at high temperatures and the formation of a Kondo resonance at lower 
temperature takes place. In addition, we find that for a dot with many overlapping resonances, an 
inverse crossover can take place. A Kondo resonance which involves many levels within the dot is 
first formed, and this coherent state is suppressed by correlation effects at lower temperatures. A 
narrower Kondo resonance, due to a single level in the dot can emerge at even lower temperatures. 



I. INTRODUCTION 



Electron-electron interactions play a crucial role in sin- 
gle electron devices (quantum dots, metallic islands), 
see for instance P, S S Hi- The electrostatic repul- 
sion between electrons in the dot can be described by 
the energy required to change the charge in the dot by 
one unit, E^- At temperatures, or frequencies, below 
this charging energy, but much greater than the spac- 
ing between the electronic levels in the dot, 6E, Coulomb 
blockade effects dominate H 0) and the conduc- 

tance through the dot is suppressed, except at degener- 
acy points. At scales below the level splitting within the 
dot, and when the ground state of the dot is degener- 
ate, electronic coherence can be restored by the Kondo 
effe ct Isl- llOlli le ading to an increase in the conductance 
[ill I12I ll.'l K^ . Hence, when the charging energy 
is much greater than the level spacing within the dot, 
Ec ^ 6E, and the number of electrons in the dot is 
odd, one expects a crossover between a high temperature 
regime where Coulomb blockade effects dominate, and a 
low temperature one described by the Kondo effect [T^ . 

The theoretical analysis of these two regimes, Coulomb 
blockade and the Kondo effect, has been carried out us- 
ing rather different techniques. Coulomb blockade can 
be studied, when the coupling between the leads and the 
dot is weak, by describing electron transport in terms of 
sequential hopping processes i-e. perturbatively in 
the tunneling matrix element. When the conductance is 
large, it is convenient to describe the internal degrees of 
freedom of the dot in terms of a collective variable, the 
phase conjugated to the total charge 0,E,llBi]. Using 
this representation, the renormalization of the effective 
charging energy by virtual fluctuations can be studied 
by a variety of methods 20,11,22, J3,J4,J5, 26). This 
approach becomes exact when the level spacing within 
the dot tends to zero and the number of transmitting 



channels is large, which applies infact to the physical sit- 
uation of metallic islands. In this limit the Kondo effect 
disappears. In the opposite limit, when the level spacing 
within the dot becomes comparable to other scales of in- 
terest, the standard procedure is to truncate the Hilbert 
space of the dot, and keep only the interacting electronic 
states closest to the Fermi level of the leads, so that the 
system is reduced to the Anderson model |0, |23|, see 
[ill (2^ (2^ . The situation when many levels within a 
magnetic impurity contribute to the Kondo effect was ini- 
tially discussed in^30]. The same effect in quantum dots 
was analyzed in[21|. Other modifications of the Kondo 
effect in a quantum dot due thepresence of many lev- 
els are discussed injUlallllllllllilllsil. Recent 
experiments!^ suggest that the regime where level spac- 
ing, level broadening and the charging energy are com- 
parable can be achieved by present day techniques. 

The aim of this paper is to analyze the different regimes 
of a quantum dot characterized by the charging energy, 
Ec, level spacing, 5E, and coupling to the leads, F, where 
F is the typical width of each level. We introduce a new 
scheme which describes the charge excitations of the dot 
in terms of a collective phase, but also retains informa- 
tion about the electronic degrees of freedom on the dot. 
This approach allows to recover the Kondo effect in the 
case of quantum impurities with exactly degenerate lev- 
els (slave rotor representation introduced in jS^I by two 
of us). Technically, the new scheme draws inspiration 
from self-consistent resummations used for describing the 
Kondo effect in single impurity models (the so-called Non 
Crossing Approximation, or NCA [iH. l4^). but adapted 
in an economical way to the situation of mesoscopic struc- 
tures (that can contain up to thousands of electrons) . 

The outline of this paper is as follows. In section ^ 
we review basic properties of the single electron box, and 
present the model that we will study throughout. The 
self-consistent approach used to analyze the model is pre- 
sented in section UllI and will be shown to describe the 
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Coulomb blockade regime of metallic grains exactly (also 
for quantum dots); we will furthermore show that it is a 
reasonable approximation for dealing with the Kondo ef- 
fect (in quantum dots). In section Hvl we will investigate 
in more detail intermediate regimes of conduction that 
were not previously attainable or considered in earlier 
works. In particular we will be interested in the different 
manner in which one goes from coherent, Kondo assisted 
transport to the Coulomb blockade regime. Extensions 
of our approach to other systems such as quantum dot 
arrays will be briefly outlined in the conclusion. 



II. MODEL FOR SINGLE ELECTRON DEVICES 

A. Phase representation and electronic degrees of 
freedom 

Let us first review the standard model for single elec- 
tron devices. The basic ingredient that determines the 



behavior of quantum dots as well as metallic islands is 



the relatively large Coulomb energy Ec 



V(2C) 



IK associated with small capacitances of the structure, 
C ~ 10~^^F. We will take also into account the internal 
structure of the small charge droplet, with energy levels 
(which are shape - and disorder - dependent), and an 
index a = 1 . . .M associated to transmitting channels to 
the outside circuit. The electrodes, to which the meso- 
scopic region is coupled, are described as Fermi Liquids 
with a spectrum Sk associated with a finite density of 
states -0(0) at the Fermi level. Physical quantities, such 
as the conductance G, are also determined as a function 
of temperature T (we set (3 = l/T) and apphed gate volt- 
age Vg. Summarizing all these energy scales, we arrive at 
the following hamiltonian: 



J 



H 



k(7 



E, 



dLd 



kap 



■ dl^aka) 



(1) 



We have introduced here the coupling between the 
mode p in the electrode and mode k in the electron box. 
The index p = 1 . . . corresponds to the energy 
levels of the dot (or metallic island), and a — l.-.M 
is the conserved channel index (including spin). We de- 
fined also Tig 8jS 0j dimensionless external gate voltage, 
Ug = eVg/{2Ec)- Because of the large number of degrees 
of freedom involved with respect to a single impurity An- 
derson model {Nl or Af can be quite large) , this hamil- 
tonian is remarkably complicated. 

One way to make progress is to single out the relevant 
charge variable on the dot, and pilot its dynamics by 
means of a rotor degree of freedom: 



^ ^ dpa-dpa 
pa 



-id/de 



(2) 



This mapping is motivated by the fact that the hamil- 
tonian -ff rotor — EcL^, with its quadratic spectrum, re- 



produces exactly the charging energy term in (Q. Al- 
ternatively, one can argue, as in [r^ll8l|. that a metallic 
grain can be studied with the same formalism used for a 
superconducting grain, except that the gap vanishes. In 
this case, the phase is the collective variable conjugate to 
the number of particles, as implied by equation 

In order to fulfill the Pauli principle, one must however 
keep fermionic degrees of freedom, so that the creation 
of a physical electron on the dot is represented as: 



dl^fle^' ^ ^ (3) 

using a channel - and momentum - carrying quasiparti- 

cle According to Q, a constraint L = J^pa fpafpo- 
must also be implemented introducing for this pur- 
pose a Lagrange multiplier S/i. Integrating explicitly the 
electrode in , one gets the following expression for the 
action of the original model in this new representation: 



J 



S= dTY,fla{dr+^p-¥)fpa + 
•^0 p„ 



^E, 



1 



13 f0 
dT / dr' 
io 



E 

pp' a 



APP\T-T')fUr)fp'a{r') e^'^(-)-^«(-') (4) 



r 



In principle, Sfj, should be integrated over [0, 27r//3] in or- One can notice in this expression how the indices are po- 
der to preserve the correct structure of the Hilbert space. 
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sitioned in the last term (cr is conserved by the charge 
transfer procedure, whereas p is not). We also have de- 
fined the bath function: 



(5) 



introducing the lead density of states D{e) = J^k ^{^"^k) 
(we have extracted here a factor l/N^ for future conve- 
nience). 



B. Metallic grains: phase only approach 



In the case of small metallic islands 0, the num- 
ber of transverse channels Af ~ 10'* — 10^ 1 is quite 
large, so that the fermionic degrees of freedom can be 
integrated out explicitely from Technically, this is 
done by expanding the initial action in the coupling 
function A^'' (r), leading to an exact expression of the 
effective action: 



S[0] = 



2 (9^0+^ 

UT 7-;^ h ingOrU 







logE^^ (n drL^ ^ A^"^'"(r™-r;„)e^'(-'")-^(-'")/,U„(r™)/,;„.„(r:,)) 

n— m— 1 ^ ^ Pm,P^<7-m 



(6) 



r 



The average over the fermionic variables ( . . . )p in the 
preceding expression is taken with respect to the free 

fermionic action Sq = Jq drj^pa fpAc>r + - S^)fpa, 
and simplifies considerably in the large M limit because 
Wick contractions at leading order occur between pair 
of fermions holding the same channel index am ■ We can 



then re-exponentiate the last term in ((BJ, and, if one 
notes that the multiplier (5/i can be treated in the grand 
canonical ensemble, i.e. reabsorbed in Ug (due to small 
charge fluctuations over many accessible states), one gets 
the effective phase-only action: 



S[0] = I 
Jo 



dr dr' ao{T—T')e 
Jo Jo 



i9(T)-ie(T') 



(7) 



r 



The kernel ao(''") in the last equation is given by the 
expression: 



ao 



pp' a 



which can be evaluated 
ao(T) 2 



at 



^sin(7rr//3)]2 
with at = NNl t^D{0)po{0) 



(8) 

(9) 
(10) 



Here po{e) — J2p^i^ ~ ^p) the grain density of states, 
which is supposed to be continous for metallic islands, 
and we have taken the large bandwidth limit to obtain 
equation 0. We have also assumed for simplicity that 
the coupling to the lead is a constant, = t (point con- 
tact), and that it is sufficiently small so that at is at 
most of order one when both Nl and A/" are large. The 
parameter at is interpreted as the (dimensionless) high 
temperature conductance between the dot and the island, 
to lowest order in the hopping t [TtI ITM . 



Physically, the use of the phase-only action I23| is vin- 
dicated by the fact that the charge fluctuations are sup- 
pressed by the Coulomb blockade phenomenon, allowing 
one to forget about the fermionic statistics of the charge 
carriers (incoherent transport through the island). 



This phase-only action O can be handled using a va- 
riety of techniques 2^, ^ ^ IHIi^ . An im- 
portant question, which is very relevant for experiments 
on small metallic grains strongly coupled to the environ- 
ment 8, 43], is the destruction of Coulomb blockade in 
the strong tunneling regime at ^ 0.1. This implies a 
generic renormalization of the Coulomb energy Ec to- 
wards smaller values, denoted here by E* the 
following, we will be more interested in this question in 
the context of quantum dots, for which the model lO is 
inadequate. However, the phase-only approach will be 
considered as an important benchmark for any approach 
that intends to describe the Coulomb blockade regime. 
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III. UNIFIED APPROACH TO SINGLE 
ELECTRON DEVICES 

A. Self-consistent scheme 

Our purpose here is to develop a simple approxima- 
tion scheme (avoiding direct numerical solution of the 
model (^) that could encompass the main interaction 
effects in single electron devices: the Coulomb blockade 
and the Kondo effect. One can note that one is usu- 
ally forced to employ different techniques to handle the 
Coulomb interaction, depending on the regime one is in- 
terested in. Since the details of the dynamics of the elec- 
trons are not usually involved in the Coulomb blockade 
phenomenon, introducing the total charge as the physi- 



J 



cally relevant variable and using the phase-only action Q 
seems mandatory to deal with this regime. On the other 
hand, dealing with the Kondo effect requires a rather 
sophisticated treatment of the correlation effects among 
electrons, and a simplification of the initial model 
is unavoidable. We will see in this paper that one can 
reconcile these two paradigms if one treats on an equal 
footing fermionic degrees of freedom and the phase vari- 
able conjugated to the total charge. 

Drawing inspiration from strong coupling approaches 
for single impurity models (NCA 4f J and its generaliza- 
tion to the phase representation |(2J) as introduced by two 
of us ^il-l^ perform a self-consistent decoupling 

of the fermion-rotor coupling in I^J: 



S ^ JdTY.fUdr + ^p-St^)fp^+ /dr|dr'^^A»'P'(r-r')(e*'^M-^(^'))/t,(r)/,,.(r') (ff) 

pp'cr 



iUgdr 



This approximation goes far beyond the lowest order per- 
turbative result that led to the phase-only action and 
is the central starting point of the present work. Its dis- 
tinctive feature is the decoupling of fermionic and phase 
degrees of freedom, whose joint dynamics is nevertheless 
determined self-consistently. Indeed, the fermionic self- 
energy obviously depends on the phase-phase correlator, 
and reciprocally. Moreover, the bosonic part of the effec- 
tive action (|I2() is similar in structure to the phase-only 
approach allowing to use the large body of work on 
this particular model. 

In order to detail the method of solution, we start here 
by simplifying the kernels appearing in the previous self- 
consistent action. The bosonic kernel can be expressed 
in terms of the full propagator of the f^^ fermions: 



= -AAi-5]A^'^'(r)Gf'(-r) 



(13) 



pp' 



where the Green's function G/ has been introduced in a 
matrix notation: 



Gf 



Gn 



Nl 



Go 



Nl 



(14) 



with the free propagator in the electron box: 



(15) 



iuj — ep + 6ii 

We also have introduced the self-energy of the quasipar- 
tides: 



(16) 



The "self-consistent phase action" that one needs to 
solve finally reads: 



dr 



+ ingdrO— d,T dr' a{T — T') 



(17) 



with q;(t) given previously. The set of equations (|13I17|) is the main technical result of this paper. In practice, 
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this is solved by an iterative procedure which starts with 
a given kernel air) as input to the action l|17|) . from 



which a new correlator ( 



') is computed. It is 



then fed back to the self-energy \i<6\ . allowing to compute 
the quasiparticle propagator (|14|l . and then a new kernel 
a(r) from equation (|13|l . This full cycle is repeated until 
convergence is reached. 

In all further calculations, we will assume a point con- 
tact between lead and box, — t, so that A^^ (r) — 

A(r) and SfV) = ^/(r) = A(T)(e'''(^)-*^(^')), which 
allows to simplify expressions and (|14|) into: 



a{T) = -J\fA{T)G'f^{-T) 



Nl 



pp' 



— E ■ 



(18) 
(19) 

^f{iuj) (20) 



Sfj, 



In the following two paragraphs, we will sketch how 
this novel scheme allows to capture both Coulomb block- 
ade and Kondo effect in single electron devices. 



B. Metallic islands: Coulomb blockade 

We consider here the case of a metallic island, for which 
the energy spectrum is taken as a continuous density of 
states /9o(e). Due to the smallness of the Fermi wave- 
length with respect to the typical transverse size of the 
island, the number of transmitting channels is usually 
quite large, Af ^ 1, as discussed before in section Hi Bl 

One interesting remark is that our self-consistent phase 
model (|llll2ll exactly reproduces the phase-only ap- 
proach (71) in this case. Indeed, we note that the 
fermionic self-energy Hlti|) is suppressed with respect to 
the kernel (|18|l by a relative factor of order 1/Af. The 
self-consistency between fermions and the phase variable 
can therefore be ignored at TV ^ 1 (one has to scale 
t cx I/VAT), leading to a free fermionic propagator: 



Gl" (r) 



itco) 



(/pV(0)/pv(r))^ 
Spp' 



(21) 
(22) 



Putting this expression back into H13|) . one recovers in- 
deed the phase-only action H7I8|) . This test case im- 
plies that our self-consistent scheme allows to deal cor- 
rectly with the Coulomb blockade phenomenon in metal- 
lic grains. 



C. Quantum dots: interaction effects 



In quantum dots, the discreteness of the energy spec- 
trum Cp and the fact that Af is generally of order one 



(unless the point contacts are quite open) invalidate the 
phase-only approach JT)). This is clear from the fact that 
coherent transport through the dot can be restored due 
to the Kondo effect at low temperature. We now sketch 
how the self-consistent action (|17|) is able to describe this 
phenomenon. 



1. Coulomb blockade regime in quantum dots 

We consider first temperatures smaller than the charg- 
ing energy Ec, but still greater than the inter level spacing 
SE. We can therefore take a continuous limit for the dot 
(due to thermal smearing of the energy levels), which 
leads crudely to the simplification: 



Mr ^ 



1 



Nl ^ ii^ - ep + 4i) 



-i7rpo(0) 



(23) 



From H20|l . this gives the 
dot: 



Gi°'=(iw) 



fermion propagator in the 



^Po(O) 



E/(icj) 



(24) 



The constant imaginary part in the last expression dom- 
inates the long time behavior of the Green's function 
G^°'^(t), so that the kernel (|18|l decays as l/r^, similarly 
to the case of the metallic island There is therefore 
Coulomb blockade in this regime, as expected. 



Z. Kondo effect in quantum dots: first discussion 

At temperatures below the interlevel spacing 5E in the 
dot, the discreteness of the spectrum becomes sizeable. 
Let us assume here that a single level sits close to the 
Fermi energy, so that we can forget all other levels. In 
this case, we can argue here that the factorization (|11I12() 
reproduces the Kondo physics correctly. This property 
is well known if, instead of the phase variable 6, one 
uses a slave boson representation of the interaction (in 
the case of the Ec — oo limit) 44]. This method has 
been extended recently by two of us to the case of finite 
Ec, using the phase representation d^^ = f^^e^^ [ioj . 
and shown to lead to a good approximation scheme for 
dealing with the Kondo effect. 

In quantum dots, the Kondo effect manifests itself by 
the buildup of a many-body resonance close to the Fermi 
level that allows the coherent transport of charge through 
the structure. The existence of the Friedel sum rule (to 
be discussed later on) guarantees that the conductance 
through the dot can recover the unitarity limit {i.e. G — 
2e^ /h) at low temperature. 

In conclusion, we can therefore expect that our self- 
consistent approach interpolates between the coherent 
Kondo regime (associated with a large and increas- 
ing conductance with decreasing temperature) and the 
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Coulomb blockade (which leads to a suppression of the 
charge fluctuations on the dot and a decreasing con- 
ductance). In practice, we still need to solve the self- 
consistent phase problem H11I12|I . This will be done in 
the next section. 



IV. INTERMEDIATE CONDUCTION REGIMES 
IN QUANTUM DOTS 

A. Solution of the phase action in the spherical 
limit 



The simplest, yet non trivial treatment of the bosonic 
action I^Tj) is to take its spherical limit ji^, @ (this 
is also equivalent to the large number of component limit 
for the sigma- model in the field theory litterature |47|). 
We introduce for this purpose the complex phase field 
X{t) EE e^^(^), with its correlator Gx{t) = {X{t)X*{0)). 
This representation is perfectly equivalent to the original 
problem if the hard constraint |X(t)P = 1 is maintained 
exactly. The approximation that we will make in or- 
der to solve the (self-consistent) phase action (|17|) is to 
impose this constraint on average only. For this we in- 
troduce a Lagrange multiplier A to enforce the equality 
Gx(-r 0) (|X(t)|2) = 1. The main drawback (on a 
qualitative level) of this approximation is that, although 
it works fine close to the center of the charge plateaus, 
it eventually breaks down at the degeneracy points [40l| . 
We will also enforce the constraint J^pa fpafpcr = L on 
average. In the symmetric case {i.e. at the center of a 
plateau), this simply implies that = 0. 

The bosonic action is now purely quadratic, leading to 
the set of self-consistent equations 



quae 



Gx{iVn) = 

Gx(r-O) = 

a(r) = 

S/(r) = 



WlT^p^/ii^n - ep) 
A{t)Gx{t) 



S/(zcj„) 



with the notation: 



A{iujn) = ^ ■ 



■ £k 



(25) 
(26) 

n -1 

(27) 

(28) 
(29) 



(30) 



We have noted here cun (resp. ;/„) a ferniionic (resp. 
bosonic) Matsubara frequency. In addition to the rela- 
tive simplicity of this scheme, we note that it allows to 
obtain a solution for real frequency quantities (doing the 
analytic continuation of these equations) , as required for 
computing spectral and transport properties of quantum 
dots. 



A few technical remarks: The system of equa- 
tions H25I3U|I was derived previously for the single-level 
Anderson model |4fl| | using a slightly different route 
(namely a large-M multichannel point of view more suit- 
able to understand some non Fermi Liquid aspects of the 
solution). We feel that a two-step procedure that starts 
with equations I|11I12|I is more appealing in the present 
context, as the self-consistent phase action (|12ll can be 
tackled in principle with variety of methods. We have 
therefore taken the point of view of simplicity in doing 
the spherical limit described in the previous paragraph. 

The numerical solution of the coupled integral equa- 
tions is straightforward using Fast Fourier Transforms 
(we will also give analytical arguments later on). In all 
calculations that follow, we take Af — 2 (single channel of 
conduction per spin). The Coulomb energy is chosen as 
the reference, Ec = 1, and the bandwidth of the continu- 
ous density of states of the electrodes is A = 50, which is 
the largest energy scale of the problem. The precise form 
of the spectrum in the electrodes playing little role in the 
low energy limit, we have chosen a semi-elliptic density 
of states: 



A{tLu) = NlIA^— [iu - iSign(c^) 



(A/2)2 (31) 



The single-level width F = |XmA(iO+)|/A^L = 4t2/A (see 
equation 10}) characterizes the strength of the coupling 
of the electrons in the dot to the reservoirs. Other im- 
portant parameters, that we will also investigate, are the 
separation between the electronic levels of the dot, 6E, 
and the number of states in the dot, N^. We also de- 
fine a typical "bandwidth" of the dot, IF = SE{Nl - 1) 
(although the spectrum is made up of dicrete states). 

A technical aspect worth mentioning is the method 
of computation of the conductance through the dot. In 
the same spirit of the decoupling performed in equa- 
tions H11I12|I . we will compute the Green's function of 
the physical electron as G^°'=(t) = G''^''{t)Gx{t). An 
analytical continuation (performed numerically) allows 
then to get the (interacting) density of states in the dot, 
p^[uj) = -(l/7r)Im G^^^iuj + iQ+). We will also use the 
interacting Landauer formula for the conductance |28ij : 



G 



Me 



-NlT / duj 



du) 



(32) 



where npiuj) = l/(e''" -I- 1) is the Fermi function. This 
expression involves the "local" density of states {i.e. 
summed over all p = 1 . . . N^) because of the form of 
the local coupling to the reservoir, see equation 

In the following, we explore situations that are inter- 
mediate to the single level Kondo effect and the Coulomb 
blockade regime, depending on the internal structure of 
the quantum dot. We insist on the fact that the results 
that we will obtain are not easily accessible to usual ap- 
proaches of the hamiltonian unless only a few levels 
in the dot are considered (for which case the Numerical 
Renormalization Group is quite successful [STlli^ ). 
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B. Results for finite interlevel spacing: Kondo 
effect and Coulomb blockade 



We now analyze the case of a few, well separated energy 
levels in the dot. We start by discussing the Coulomb 
blockade using the spherical limit (for temperature larger 
than the interlevel spacing), then we give analytical ar- 
guments in favor of the existence of a Kondo resonance at 
lower temperature, and finally wc show the full numerical 
solution of the self-consistent equations corresponding to 
the usual situation found in experiments. 



1. Renormalized Coulomb energy 



As discussed in section ITlI C II for T ^ SE the quan- 
tum dot is in the Coulomb blockade regime, and the 
bosonic kernel H28|l behaves at long times as a(r) ~ 
at/r^, see eq.Q. Here at is a measure of the dimension- 
less high temperature conductance. One can now focus 
on equation (|26|l . that we write (at low temperature) as: 



1 

2^ t/2 /{2Ec) + Trat\u\ + X- a{iO) 



1 



(33) 



The renormalized Coulomb energy E* is obtained as the 
mass term in the phase propagator H25|) . which decays 
exponentially over time scales of order Ti/E*. Solving 
the previous equation, one gets: 



e: = 



A-a(iO) 



Trat 



27rat Ecc-^ "* 



(34) 



which is exponentially small in the bare conductance at- 
This quick calculation allows us to understand the ori- 
gin of the Coulomb blockade for quantum dots in our 
formalism. 



2. Single level Kondo effect: analytical proof 

In quantum dots, the basic manifestations of the 
Kondo effect are twofold: existence of a small energy 
scale, the Kondo temperature Tjf, and restoration of 
the unitarity limit (at T < T^), as discussed in sec- 
tion UTEni 

To prove the first point from our integral equations \25l 
I30|l . we examine qualitatively the solution of the self- 
consistent equations at temperatures smaller than the 
interlevel spacing 6E. There, the fermionic Green's func- 
tion H27|l is dominated by the ep — pole, so that 
G^f^ir) - ~1/{2Nl) Sign(r), at long times. From 
equation if^ . this produces in turn a rotor self-energy 
Q!(w) -(27Vr/7r) log li^/Tol at low frequency (here To 
is an undetermined cutoff, in practice of order ^/Ec^). 
Upon lowering the temperature, this self-energy can 
reach the charging energy E^. (assumed to be only weakly 
renormalized), indicating the suppression of Coulomb 



blockade and the beginning of the Kondo regime. This 
happens (very roughly) for -(27Vr/7r) log(TK/To) ~ E^, 
giving the estimate: 



TK{NL)^Toexp - 



ttE, 



2NT I 



Tociq,|-|^~| (35) 



This is the well-known value of Tk in the local moment 
regime (at half- filling) , if one uses the standard notation 
Ec = U/2 with Af = 2 (single-channel). It clearly ex- 
plains that one key step to the experimental observation 
of the Kondo effect in quantum dots Ts] lies in the real- 
ization of strongly coupled structures, having T compara- 
ble to Ec in magnitude, such that the Kondo temperature 
remains accessible. 

We finish by explaining precisely the restoration of full 
coherence below the Kondo temperature. In order to 
do this, we set the temperature to zero and perform an 
exact low-energy solution of the system of equations H25t - 
I3U|I (this was done in Appendix-C of ^3)- This analysis, 
valid while SE for an arbitrary number of levels, 
leads to the following value of the zero temperature, zero 
frequency density of states in our approximation: 



Pdiuj = T = 0) 



1 



tt/2 



ttNlT TV -t- 1 



tan 



(^ 



7r/2 
JV+1 



(36) 



The interpretation of this relation is that, whenever the 
temperature is lower than Tx, the density of states is 
pinned at its non-interacting value, no matter how large 
the Coulomb energy Ec is. This reflects the presence of 
the Kondo resonance at low energy. The exact Friedel's 
sum rule however is p^^'^^{uj = T = 0) = l/(7rA^Lr), in- 
dependently of Af. The difference between equation H36(l 
and this exact result is a consequence of the decoupling 
approximation made above, and is related to the non 
Fermi Liquid features described in ref. ^^l- This arte- 
fact is nevertheless quite small in practice, since goes 
to the exact value for Af large, and is only 10% off for 
Af = 2. The Friedel's sum rule provides also a simple 
explanation for the restoration of the unitary limit in 
the conductance below the Kondo temperature, since the 
Landauer formula l|32|) at zero temperature leads to: 



G(T = 0) = 



Afe 



■[T:NLTpd{Q)\ 



(37) 



3. Kondo effect: numerical solution and crossover to the 
Coulomb blockade regime 

The numerical solution of the set of self-consistent 
equations H25I3U|I allows in principle to investigate all 
regimes of parameters. We will first concentrate on the 
gradual suppression of the Kondo effect with decreasing 
values of the single-level coupling F. In order to maintain 
the Coulomb blockade effect, we will keep in this section 
the total (multilevel) coupling fixed, F"^'" = iV^F = 0.5, 
and vary the number of levels to allow changes in 
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r = r'""i*7A^L. We will also fix the total bandwidth of 
the dot, W = 1, so that the level spacing (assumed to 
be uniform) is also decreasing, 6E = W/{Nl — !)• This 
way of proceeding allows to interpolate from the few, 
well-sepatated levels situation relevant for small quan- 
tum dots in the Kondo regime to the case of larger dots 
with small level spacing and many levels N^, which shows 
only Coulomb blockade. 

The figurendemonstrates indeed how the low temper- 
ature local density of states evolves from Nl = 1 (single 
level: regular Kondo effect) to Nl = oo (continuum of 
levels: Coulomb blockade only). In particular the rapid 
suppression of the Kondo peak for diminishing values of 
r at increasing Nl is in accordance with our previous 
discussion of the Kondo temperature, equation H35|) . 




FIG. 1: Local density of states pd(w) at Ec = 1, (5 = 800, and 
fixed W = {Nl- 1)SE = 1 and T"""'*' = NlF = 0.5. The four 
curves correspond to different values of the number of levels 
in the dot: Nl = 1,3, 5, oo (from top to bottom, following 
the evolution of the central peak). 

The temperature dependence of the electronic spec- 
trum is presented for the three level case Nl = 3 in fig- 
ure |21 When temperature is lowered, the zero frequency 
density of states starts diminishing (by Coulomb block- 
ade of states with different charge). One then reaches a 
minimum, before Pd{0) begins shooting up, towards the 
unitarity limit (Friedel's sum rule) at zero temperature. 

It is useful to compare this evolution of the density of 
states to the variations of the conductance G(T, Nl) with 
temperatures, figure |3 This figure illustrates the reduc- 
tion of the Kondo temperature with T by the downward 
shifting of the minimum of conductance. The Coulomb 
blockade is present at higher temperature, as shown by 
the decrease of G{T) for Tk < T < E* upon lowering T 
[TgI I . For the last curve with a continuum of states in the 
dot, this behavior persists up to zero temperature. The 
inset in log-scale on the same plot allows to grasp more 
clearly the saturation of conductance below the Kondo 
temperature. 
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FIG. 2: Temperature dependence of pd(t^) for the case Nl = 
3; temperatures correspond to /? = 800 (sharp central peak), 
/3 = 25 (broad peak), /9 = 10 (deep minima) and P = 2 
(shallow minima filled by thermal excitations). 
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FIG. 3: Conductance G{T, Nl) in units of /h for the same 
parameters as in figure Q in function of temperature; curves 
with Nl = 1,3,5,00 follow from top to bottom in the ex- 
treme left of the plot. Inset: G{T,Nl = 1) in temperature- 
logarithmic scale. 



C. Effects due to overlapping resonances in 
multilevel dots 

The analysis in the previous section describes the usual 
situation where a crossover from the Coulomb blockade 
regime to the Kondo effect takes place as the tempera- 
ture is lowered below the interlevel spacing in the dot. 
A different behavior can be expected when there is a set 
of overlapping resonances at low energies within the dot, 
i.e. an ensemble of levels of individual width greater than 
their separation, T ^ 6E, that act together as a single ef- 
fective level with enhanced coupling to the leads (the typ- 
ical bandwidth of this set of level should also be smaller 
than the charging energy). The presence of broad reso- 
nances near the Fermi level can be relevant to some ex- 
perimental situations 35, '36^. Note that the conductance 
distributions for open ballistic quantum dots, estimated 
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from Random Matrix Theory are quite broad ^^IH^I- We 
give first some qualitative arguments in order to describe 
the new effects expected in this regime, and then present 
explicit calculations using the integral equations. 



1. Inverse crossover from the Kondo effect to Coulomb 
blockade 

In practical situations encountered in quantum dots, 
the single-level Kondo temperature Tk is usually much 
smaller than the level spacing 5E (which is typically com- 
parable to Ec). Therefore Coulomb blockade occurs (if it 
does) inevitably before the Kondo effect sets in, as was 
shown at length in the previous section. There is however 
a simple mechanism that allows to obtain Kondo temper- 
atures greater than both level spacing and renormalized 
Coulomb energy E* . The idea is that when the individual 
level width T exceeds the inter level spacing 5E^ many lev- 
els are involved in the formation of the Kondo resonance. 
This leads to a "multi- level Kondo temperature" [s^, |^ : 

which can be greatly enhanced with respect to the single 
level estimate, equation (|35|) . by the presence of many 
levels Nf.ff > 1 acting together (increasing the number 
of channels Af might also contributes to this effect). This 
way of enhancing the Kondo temperature allows to ob- 
tain a new regime where Tj^"'*' ^ E*,SE, so that the 
Kondo effect can now occur before the Coulomb block- 
ade (when lowering the temperature), in an inverse man- 
ner as observed traditionally in quantum dots. The fact 
that the Coulomb energy can be strongly renormalized 
to smaller values adds credibility to this idea. One fur- 
ther notes that, at even lower temperatures, a Kondo 
peak associated to the formation of a resonance which in- 
volves only one electronic state in the dot will ultimately 
emerge. One therefore has a "two-stage Kondo effect" 
(if the single-level Kondo temperature is not vanishingly 
small) . 

In order to be more precise, we will first give a concrete 
example with a limiting case that one can understand 
independently of any approximation scheme. Then, we 
will illustrate in detail this "inverse crossover" using our 
integral equations. 



2. Limit of exactly degenerate levels 

We consider here the extreme limit in which Nl levels 
in the dot are simultaneously put to zero: Cp = for all 
p (the total bandwidth W is therefore also zero) . We can 
formulate the model after a redefinition of the fermionic 
operators (unitary transformation) {dpg.} ~* {'^po-}' such 



that: 

^L-7^E< (39) 

In this case, the remaining fermionic degrees of free- 
dom, cj,^ for p > 1, simply decouple from the problem, 
leaving an effective single level Anderson model describ- 
ing the fermion c\^. Actually, a capacitive coupling per- 
sists between this fermion and all the other ones, but this 
influence gets frozen at low temperature. As the con- 
ductance through the dot is obtained from the Green's 
function: 

Grir) EE i-^(4,(0)dp,.(r)) = (cL(O)ci.(r)) 
pp' 

(40) 

one gets a full Kondo effect and a complete restoration of 
unitary conductance at low temperature. Furthermore, 
the width of this effective level is simply iV^r, as one 
checks by inserting in equation This leads then 
to an enhanced Kondo temperature, 

as discussed in the introductory part of this section (be- 
cause 6E = in this limiting case, one has A^e// — Nj^). 

We can easily check that our self-consistent scheme 
preserves this interesting property of the model. Indeed, 
when all levels €p are exactly degenerate {W = 0), one 
gets from equation (|27|l the /-electrons Green's func- 
tion Gf{iuj) = [ioj — I]/(«a;)]~^. We obtain therefore the 
Kondo effect of a single effective level 40]. 

The following paragraph will allow us to make this dis- 
cussion more meaningful, by studying the more realistic 
case of a quantum dot in the regime T ^ SE, correspond- 
ing to nearly degenerate levels. 

3. Inverse crossover: illustration 

We now illustrate the inverse crossover discussed quali- 
tatively in paragraph llV (] Ti bv solving our integral equa- 
tions in the regime F ^ (5E, where multi-level effects play 
an important role. We will assume here that T <^ Ec 
so that one can neglect the single level Kondo effect at 
low temperature (see however the next paragraph). For 
this computation, we have fixed F = 0.04, and taken 
Nl = 9 states in the dot, varying the interlevel spacing 
from 6E — (exactly degenerate level case considered in 
the previous paragraph) to SE = 0.01 ^ F (possibility of 
multilevel effect) to SE — 0.08 > F (absence of multilevel 
effect). 

The low-temperature local density of states displayed 
in figure 01 shows the expected multi-level enhanced 
Kondo peak a,t 5E = corresponding to formula l|41|l . 
Upon increasing the level spacing to SE — 0.01, Coulomb 
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blockade sets in at a scale E* < 



herence effects remain around T]^"'* 



rpmulti 



however co- 
(since we are in a 
regime with 5E < T). This results in a surprising split- 
ting of the Kondo resonance at low energy. The last 
curve is taken with 5E — 0.08 > F, so that no multilevel 
Kondo effect is possible, and only Coulomb blockade is 
observed. Another interesting consequence of this phe- 



0.8 r 



havior of the /-Green's function at high temperature: 




FIG. 4: Density of states pdi'jj) on the dot for Ec — 1, 
F — 0.04, inverse temperature (3 = 400, Nl ~ 9 states, and 
different values of the interlevel spacing: 5E — 0, 0.01, 0.08. 

nomenon is that the temperature dependance of the con- 
ductance is reversed with respect to the usual signature 
of the Kondo effect in quantum dots, i.e. to figure O 
Indeed, upon lowering the temperature, one notices an 
initial increase of the conductance (due to the multilevel 
Kondo effect), then a sharp decrease of the conductance 
because of the Coulomb blockade. This is illustrated by 
the middle curve in figure |31 




G 



(42) 



Equation then leads to |58l |: 
NT 



51 log- 



(43) 



The Kondo temperature is reached when this kernel is 
of the order of Ec^ so that one finds the final equation 
which determines T™'*': 



n 



K 



(rr"')2 + (kp| 



To 



(44) 



where Tk is the single level Kondo temperature H35|l . 
A similar result was obtained previously by a renormal- 
ization group argument 'SOj. The limiting cases stud- 
ied before are obviously contained in the previous equa- 
tion: T'jJ'"'*' reduces to Tk for widely separated levels 
{6E large) and in the opposite limit of exactly degener- 
ate levels (or ii W <^ T™""'), ^1?""' = To{Tk /ToY^'^^ 
consistently with equation H41|l . In general, T™^'*' is en- 
hanced with respect to the single level Kondo temper- 
ature Tk. As an example, formula 1441) shows that in 
the regime SE <C Tg"''*' < TF^ < iJc, one obtains then 
jnmuiti ^ To exp[-7r£:c<5E/(2A/'Frf "'")], so that the num- 
ber of effective levels taking part to the multilevel Kondo 
resonance is Neff = T^^^^^/SE, as discussed qualitatively 
in equation 

We conclude this paragraph by summing up the physi- 
cal picture that leads to the observed non-monotonous 
conductance. In the case of many overlapping reso- 
nances, a quantum dot can be described as a small metal- 
lic grain dominated by Coulomb blockade at low temper- 
ature, which implies a vanishing zero-frequency density 
of states. Upon raising the temperature, many different 
charge states become available by the thermal smearing 
of the Coulomb blockade, and the conductance is rapidly 
increasing on a scale of the order of E* (which is also 
the typical size of the dip observed in the split Kondo 
peak). Due to the large single level width considered 
in this regime, all these energy levels can then act coher- 
ently as a localized spin degree of freedom that undergoes 
the Kondo effect. This explains the upturn of the con- 
ductance when temperature reaches the Kondo energy 

T^multi 



FIG. 5: Conductance G{T) for the same parameters as in 
figure 01 

One can also perform a general evaluation of the mul- 
tilevel Kondo temperature, starting from the leading be- 



4. Two-stage Kondo effect 

Finally we consider again a multilevel case, (SE <C F, 
but now with T ^ Ec, so that the single-level Kondo reso- 
nance is accessible to the low temperature regime. There- 
fore, one witnesses a further increase of the conductance 
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at low temperature, taking place after the inverse Kondo- 
to-Coulomb crossover that we discussed previously. The 
occurence of such a "two-stage Kondo effect" is depicted 
for the conductances shown in figure El For this calcula- 




'o 0.1 0.2 0.3 

T 



FIG. 6: Conductance G(T) in the case Nl = 5, T ^ 0.1 and 
various level spacings: 5E = 0, 0.025, 0.05, 0.1, 0.25. 

tion, we have taken N]^ — 5 levels, F = 0.1 and various 
level spacings between zero and 2.5F. 

The two curves with < SE < T show indeed this 
two-stage Kondo effect: a first rise of the conductance 
at high temperature due to the multilevel resonance at 
T - Tg^'"'**, then the Coulomb blockade at T - E* and 
then a further increase at T ~ Tk (smallest scale of the 
problem). The last two curves, with 6E = 0.1, 0.25, show 
the usual crossover that was illustrated on the figure 13 in 
section HVBl 

A last point that was checked is that our results arc 
weakly sensitive to the addition of more levels outside 
the energy window of width F. This calculation is shown 
in figure [7| 



5. Summary of the different regimes of transport 

A sketch of the different regimes analyzed in this paper 
is given in figure|Hl as a function of the level spacing, SE, 
and single level width, F. We assume that the number 
of electrons in the dot is odd, so that the ground state, 
in the absence of the coupling to the leads is degener- 
ate. The transition between different regimes is a smooth 
crossover. We will not discuss in the following the effect 
of the renormalization of the charging energy, therefore 
the zone boundaries that locate the different regimes will 
not be determined here. However, this can be done in 
practice by solving the integral equations, as discussed 
in the body of the paper. Also it would be interesting 
to perform a full Renormalization Group analysis of the 
renormalization of both and E* along the lines of 

Pill 
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FIG. 7: Lower curve: conductance G{T) in the case ~ 5, 
r = 0.1 and SE — 0.015. Upper curve: similar model, but 
10 additional levels (with larger spacing SE' — 0.2) have been 
superimposed to the previous ones. 

The occurence of Kondo effect is signaled by a non 
monotonous temperature dependance of the conduc- 
tance, and is observable at F <C for temperatures 
smaller than the single-level Kondo temperature Tk 
given by equation (|35|l . Depending on the relative val- 
ues of level spacing SE and Coulomb energy Ec, Coulomb 
blockade might also be present above the Kondo temper- 
ature, which is the usual experimental situation encoun- 
tered in quantum dots. We emphasize that if F is much 
smaller than E^. the Kondo temperature is extremely 
small and the Kondo effect inobservable in practice; this 
corresponds to weakly coupled dots, and this situation 
shows only Coulomb blockade. In the case SE T we 
have defined a region "MultiLevel" which corresponds ac- 
tually to various regimes discussed previously. This de- 
scribes the inverse crossover (Kondo to Coulomb Block- 
ade), as well as the two-stage Kondo effect shown in sec- 
tion IIV (T4l This region can also imply that the Kondo 
effect is observed in the usual manner, but with a Kondo 
temperature greater than the single-level estimate. Note 
that the Kondo effect involving many levels can also oc- 
cur in dots with an even number of electrons, in a sim- 
ilar manner to the Kondo resonance which arises at a 
singlet-triplet crossing in an applied magnetic field [s^. 
Finally, the region F ^ of the diagram (denoted "Co- 
herent") is associated with large conductances that are 
weakly modulated with temperature or applied gate volt- 
age. 



V. CONCLUSIONS 

In this paper we have presented a method of calcula- 
tion for strongly correlated mesoscopic systems in terms 
of a collective phase variable and the quasiparticle de- 
grees of freedom. The scheme is valid both for the study 
of the Coulomb blockade regime and of the Kondo effect. 
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1 hElEc 

FIG. 8: Sketch of the the different regimes discussed in the 
paper, as function of the charging energy, Ec, level spacing 
within the dot, 5E, and level widths, F, with a constant num- 
ber of levels. The notation "Kondo" corresponds to a single- 
level Kondo effect. The region "Kondo -|- Coulomb-Blockade" 
is the usual situation in quantum dots. The "Multilevel" re- 
gion is also associated to the Kondo effect, but with impor- 
tant renormalization of the Kondo temperature, as discussed 
in the text. The "Coherent" regime stands for a temperature 
independant conductance. 

We have shown examples of the crossover between 
a Coulomb blockade regime at temperatures below the 
charging energy, and the formation of a Kondo resonance 
at temperatures lower than the separation between levels 
within the dot. In addition, we have described an inverse 
regime, where a Kondo like resonance is splitted at low 
energies by Coulomb blockade effects. This regime is as- 



sociated to the existence of many conduction channels, or 
overlapping resonances within the dot, which contribute 
collectively to the Kondo effect. The coherence of this 
state is destroyed at low energies by Coulomb effects. At 
even lower energies, a narrow Kondo peak, asociated to 
a single level within the dot, will emerge. 

Because the self-consistent approach used here was suc- 
cesfully applied to model of strongly correlated electrons 
(Hubbard model) in a previous work ^ we can also en- 
vision possible applications of this work to the physics of 
granular materials or quantum dot arrays in the vicinity 
of the metal-insulator transition [5a . l54Ll55ll5q| . Disorder 
effects, that were neglected here, should also be included 
in future work along these lines. 

Other possible development is an exact (numerical) 
solution of the self-consistent action H11I12|I instead of 
taking the spherical limit as was done in section IIV Al 
which brings some limitations (as the restriction to val- 
leys of conduction). Also importantly, we believe that 
the present work, both by the technical concepts and the 
physical ideas, connects in an original manner the fields 
of strongly correlated systems and mesoscopic physics. 
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